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Abstract 

We present, in the framework of the canonical quantization, a class of nonlinear Schrodinger 
equations with a complex nonlinearity describing, in the mean field approximation, systems of 
collectively interacting particles. The quantum evolution equation is obtained starting from the 
study of a iV-body classical system where the underlined nonlinear kinetics is governed by a kinetic 
interaction principle (KIP) recently proposed [G. Kaniadakis: Physica A 296 (2001), 405-425]. 
The KIP, both imposes the form of the generalized entropy associated to the classical system, 
and determines the Fokker-Planck equation describing the kinetic evolution of the system towards 
equilibrium. 
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The Schrodinger equation is one of the most studied equations both from a mathematical 
and a physical point of view. A particular interest concerns the possible nonlinear exten- 
sion of this equation. Just one year after its discovery, Fermi proposed the first nonlinear 
generalization 

Among the many nonlinear extensions of the Schrodinger equation (NSE) it is worthy 
to recall the effort made by Bialynicki-Birula and Mycielski (BBM) ^2 1 who proposed to 
introduce the simple nonlinearity —b ln(|^| 2 ) ip, selected by assuming the factorization of 
the wave function for composed systems. Motivated by this work, some experimental tests 
were carried out based on neutrons interferometry 3J], and Fresnel diffraction with slow 
neutrons |4|. They gave the upper limits b < 3.4 x 10~ 13 and b < 3.3 x 10~ 15 , respectively, 
suggesting that there is no a real ground for a such nonlinear term in the Schrodinger 
equation. 

Weinberg in order to save partially the superposition principle of the standard quantum 
mechanics, proposed a very general class of NSEs with a homogeneous nonlinearity. A 
further experimental test p] shows again a sharper bound on the deviation from linearity. 

Notwithstanding, many complex systems can be described by introducing nonlinear terms 
in the linear Schrodinger equation, to take into account the interactions among many par- 
ticles. Applications of the BBM equation to many particles system, in particular in the 
nuclear physics, have been considered in |7( (and references therein). There, it has been 
suggested that, in spite of the negative tests based on neutrons, the BBM equation and 
more generally, NSEs, can be usefully applied in the description of extended objects. 
The NSE with the cubic nonlinearity [a, |i| has been used to study the dynamical evolution 
of a boson gas with a 5-function pair- wise repulsion or attraction |l0| , and in the description 
of the Bose-Einstein condensation of alcali atoms like 7 Li, 23 Na and 87 Rb [llj. In |l2j a 
NSE with complex nonlinearity has been introduced in order to account for a generalized 
Pauli exclusion- inclusion principle between the quantum particles of the system. Finally 
it is worthy to mention the family of the Doebner-Goldin equations Q], introduced from 
topological considerations, as the most general class of Schrodinger equations compatible 
with the Fokker-Planck equation, for the probability density p = \ip\ 2 - 

The BBM equation is closely connected with the Boltzmann-Gibbs (BG) entropy. In 
fact, it can be shown that, within this theory, the binding energy released in the splitting 
of an arbitrary wave function, for an isolated system, in non-overlapping parts of the same 
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form ip(x) — > Y,iP]^ z '">p{.x — x.), isgiven by AE = — b Y^iPi m Pi? where p i are normalized 
according to J2iPi — 1- I n R- e f- ij it has been suggested that the BBM equation can be 
used to formulate the thermodynamical description of spatially extended or self-interacting 
quantum mechanical objects. 

In the present paper we approach to the canonical quantization of a classical system, 
where the statistical information is supplied by a very general entropy. 

We recall that, in the recent years nonlinear generalizations of the Fokker-Planck equation 
and their connection with generalized entropies have attracted the interest of many authors 

bremsstrahlung 3], self-gravitating systems [3| 
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few examples, among many others, where we can find experimental evidence for the relevance 
of the use of generalized entropies and nonlinear Fokker-Planck equations (NFPEs). 
In Ref. Q, 

it has been proposed a kinetic interaction principle (KIP) which defines a 
special collective interaction among the iV-identical particles of a classical system. The KIP 
imposes the form of the generalized entropy and governs the evolution of the system toward 
the equilibrium, by fixing the expression of the nonlinear current appearing in the NFPE. 

We present the Lagrangian formulation of the model. By using the minimum action 
principle, we obtain a class of nonlinear Schrodinger equations containing a complex 
nonlinearity and describing, in the mean field approximation, a quantum system of 
interacting particles. The form of the entropy of the ancestor classical system determines 
the nonlinearity in the evolution equation. Finally, by means of a gauge transformation 



of third kind 



13, 



22| we show that this family of NSEs is gauge equivalent to another 



family of NSEs having a pure real nonlinearity depending only on the probability density 
field. 

Some evolution equations obtained from kinetic equations known in literature within 
the Boltzmann-Gibbs entropy, the Tsallis entropy, the quantum-group entropy and the 
Kaniadakis entropy are given to illustrate the relevance of the method. 



Let us consider a classical Markovian process in a n-dimensional space 
dp{t,x) 



dt 



ir(t, y -> x) - n(t, x -> y) dy , 



(1) 



describing the kinetics of a system of N identical particles. We assume for the transition 
probability 7r(t, x — > y) a suitable expression in terms of the populations of the initial site 



x and the final site y. 

According to KIP [20] we pose 

ii(t, x -> y) = r(t, x, x-y) a(p) b(p') c(p, p) , (2) 

where p = p(t, x) and p' = p(t, y) are the particle density functions in the starting site x 
and in the arrival site y, respectively, whereas r(t, x, x — y) is the transition rate which 
depends only on the starting x and arrival y sites during the particle transition x — > y. 
The two factors a(x) and b(x) in Eq. (J2J) are arbitrary non-negative functions, supposed to 
be continuous. Their limit values for x — > are subject to precise physical requirements. 
For instance, a(0) = because if the starting site is empty the transition probability is equal 
to zero; b(0) = 1 because no inhibition occurs if the arrival site is empty. Finally, the last 
factor c(p, p') takes into account that the populations of the two sites can eventually affect 
the transition, collectively and symmetrically. 

The expression of the functional b(p) plays a very important role in the particle kinetics 
because can stimulate or inhibit the transition x — > y allowing in this way to take into 
account interactions originated from collective effects. 

With the assumption given in Eq. (J2J) for the transition probability, under the f irst 
neighbor approximation, we can expand the r.h.s of Eq. (0) up to the second order in y 



dp d 
dt dx. 



'Q+ 9 S 7(p) + C„7(p)g|ln K (p) 



(3) 



with % = 1, ■ ■ ■ , n and summation over repeated indices is assumed. 

In Eq. (JHJ) 7(p) = a(p) b(p) c(p, p) and k(p) = a(p)/b(p), whereas the coefficients Q and Q, 
are given by 

C = / x , y) Vi d v 5 and C y = g / x > y) Vi Vi d y ■ ( 4 ) 

In the following we consider the case of linear drift by posing j(p) = p and after intro- 
ducing the quantity thrift L = — C< — ^Q./dx., the z-th component of the drift velocity and 
assuming £ y = D8 ij} with D the diffusion coefficient, Eq. (JHJ) becomes 

+ V • {p [n drift - DV ln«:(p)] } = , (5) 

which is a NFPE for the field p. The total current j = j'drift+idiff i s the sum of two quantities: 
i drift = P '"drift which is the standard linear drift current and j dm = — Dp Vln/t(p) which is 
a nonlinear diffusion current different from the standard Fick's one. 



4 



Eq. (jSJ) describes a class of nonlinear diffusive processes varying the functional n(p). In 
the following we assume a constant diffusion coefficient D = const. 

The functional n(p) is related to the entropy of the classical system as can be showed in 
the non-equilibrium thermodynamics framework 2, 24], by relating the production of the 
entropy to a NFPE. In fact, we may show this link on a more general ground by starting 
from a constrained entropic form 



S(p) = S(p)-(3 J S(x)pdx-(3' J pdx , 



(6) 



where (3 and (3' are the Lagrange 
E = J S(x) pdx and the normalization J pdx = 1. 
Without loss of generality, we can assume 



multipliers associated to the constraints of the total energy 
■ , : - ' - ix = 1. 

for the entropy function the following expression 

S(p) = -k J dx J dp ln«(p) . (7) 



The standard BG-entropy is recovered from Eq. (J7J) after posing n(p) = e p. 

Quite generally, the evolution of the field p in the configuration space is governed by the 
continuity equation 



d p 
~dt 



+ V- J = 



(8) 



and we assume a linear relation between the current J and the constrained thermodynamic 
force V(5S/5p) 



T D t SS 



Putting Eq. Q into Eq. (jSJ) we obtain 



dp «-, i 



— (3V £(x) - DV Innip) 
k 



and, after introducing the drift velocity 



D 



"drift = (3V £(x) 



(9) 



(10) 



(11) 



Eq. (jTOjl takes the same form of the continuity equation (0). 

Let us now introduce the wave function if> describing, in the mean field approximation, 
the quantum analogous of the classical iV-particle interacting system depicted above. We 
postulate that the evolution equation of this system is given by a NSE 



i h 



dip 
~dt 



2 m 



Aip + A(ip*, ip)ip + V(x) ip 



(12) 



where A(if)*, if)) is a complex nonlinearity and V(x) is the external potential. 

In the canonical picture, the motion equation (j!2|) can be obtained after introducing the 
action of the system 

A= f C(ip\ ip)dxdt , (13) 



where £ is the Lagrangian density given by 



vo = i \ U* ^ - ^ A - ^ i v i>\ 2 - u{p, ip) - r v( x ) ip , (i4) 



and U(if)*, if)) is a nonlinear potential describing the collective interaction between the par- 
ticles. 

The NSE ()12j) is obtained by applying the minimum action principle to the system de- 
scribed by the Lagrangian (fT4|) 

14 = 0. (15) 

5 if)* v ' 

The nonlinear functional A (?/>*, if)) in Eq. (fT2*|) is given by 

HP, V0 = t4t [U{r, pdxdt. (16) 



Sip* 

It is convenient at this point to introduce the hydro dynamic fields p and E related to the 
^-function through the polar decomposition 25|, |26j 

if(t, x) = p x '\t, x) exp (i E(t, x)) . (17) 

Introducing the real part ?/>) and the imaginary part W(ip*, ip) of the complex non- 

linearity A(ip*, ip) = W(ip*, ip) +iW(ip*, ip), Eq. (JT2J1 can be separated in a coupled of real 
nonlinear evolution equations for the phase E 



and the density p 



|£ + V.(^p)-|p WteS) =0. (19) 



In Eq. fll9j) j = (V E/m) p is the linear current of the standard quantum mechanics. 
Both Eq.s (fT%|) and (fTT)|) can be obtained from the variational problem 

^ = 0, and 6 -^ = 0, (20) 
op oh 
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where the Lagrangian density (JH|, after taking into account Eq. (JT7J), assumes the expres- 
sion 



C(p, E) 



<9E (VE) 2 ft 2 (Vp) 2 



f/( P , s)-y(*) P . 



(21) 



dt r 2m »m p 

We observe that the nonlinear functionals W(p, E) and W(p, E) are related to the nonlinear 
potential U(p, E) through the relations 



S 



h 5 



W(p, E) = — / U(p, E) <fe dt , W(p, E) = — / [7(p, E) dx dt , (22) 
dpi 2 p o E 7 



and by choosing [7 (p, E) as 



U(p, E) = — D p V ln«(p) ■ VE 



(23) 



from Eq.s (|18|) and (|19|) we obtain 



<9p 

~dt 



+ V 



VE 



m 



p — D p V In k(p) 



. 



(24) 
(25) 



By posing thrift = V E/m, the quantum drift velocity, Eq. (|23j) assumes the same expression 
of the Fokker-Planck equation © describing the kinetics of the classical system obeying to 
KIP. 

The NSE in the ^-representation is given by 



dtp 



h 



d 



%D 



Aip + Dp — In /<p)AEV + i 
2m ap 2p 



p Vlnft(p) V + V(a;) ^ 



(26) 



where the nonlinearities W and W are determined through the functional n(p) which also 
defines the entropy (0) of the ancestor classical system. 

nn 

We introduce now the following nonlinear gauge transformation of third kind |21|, |22J 

(27) 



ip — ► (p — ip exp f — — m D In «(p) 

which, being an unitary transformation, does not change the amplitude of the wave function, 
l^| 2 = I^P = Pi an d transforms the phase E of the old field ip in the phases a of the new 
field according to the equation 



a = E — m D In «(/?) . 



(28) 
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After performing the transformation (|27p . Eq. ()26|) becomes 



d(j) ft 2 id 1 r -|2~ 

i h —— = A + m D 2 < p — — In /t(p) A In nip) H — V In 

at 2m I ap 2 L J 



0+^(^)0 



(29) 



where in Eq. f!29|) appears a purely real nonlinearity depending only on the field p. 



To illustrate the relevance of the approach described above, we derive some NSEs starting 
from as much entropies known in literature: BG-entropy, Tsallis-entropy, quantum-group 
entropy and K-entropy. 

(a) As a first example we consider the standard classical entropy 



£bg(p) = —k J p\npdx 



(30) 



which follows from Eq. (J7J) by posing a(p) = e p, b(p) = 1 so that n(p) = e p. The 
corresponding continuity equation becomes 



that is the linear Fokker-Planck equation ^J. 

The NSE is readily obtained from Eq. (|29jl and becomes 

-2 



(31) 



at 



2m 



Acb + mD 2 



Ap 




!?)*] 


P 







which was studied previously in 
Schrodinger equation 



27| . In particular, Eq. 



dt 



r 



2m 



A X + V(x) X 



<P + V{x)<f), (32) 
is equivalent to the linear 

(33) 



with % = h J 1 — (2mD/h) 2 and the field x 1S related to the hydrodynamic fields p and a 
through the relation x — P 1 ^ 2 exp(icr/£). 



(b) As a second example we consider the Tsallis entropy 



S q (p) = k I ^—^-dx 



l-q 



Eq. (J34j) can be obtained from Eq. (J7|) after posing 



(34) 



In k(p) = — ln 9 



a 



(35) 



where a = gV(i-<?) is a constant and the deformed g-logarithm is defined by 

hi fx) : 



1-q 



(36) 



which reduces to the standard logarithm for q — > 1 (ln^x) = lnx) as well as Eq. 
reduces to the BG-entropy in the same limit. 

The evolution equation for the field p is given by 



at \ m 

which is the NFPE discussed in Ref. HE]. 

The NSE associated to Eq. (J37)) assumes the form 



, 



(37) 



i h 



~dt 



h 2 



2 m 



A(f> + q 2 mD 2 p 2 ^ 



Ap 
P 



4> + V(x)(f), (38) 



and reduces to Eq. (j32jl for q — > 1. 

We observe that the nonlinearity in Eq. ()38)1 coincides with the same one appearing 
in the NSE proposed in Ref. j^J, after replacing q — ► 2 — q, and obtained in the Tsallis 
statistics framework by using a different approach. 

(c) Another interesting example is provided by the quantum groups entropy (jyj 



S q (p) = -k [ pQ P \ dx 

q J Q — Q 1 



obtained from Eq. Q with the choice 



In n(p) 



qpo- 1 - q- 1 p q 



(39) 



q-q- 



(40) 



Entropy (}39|) is invariant under the interchange q «-> q^ 1 and reduces to the BG-entropy in 
the q — > 1 limit. 

The evolution equation for the field p is given by 



dp 
~dt 



+ V 



VE 



rn 



p-DV 



q - P « 



q + l 







(41) 



and the associated NSE is readily obtained from Eq. ()29j) and assumes the form 



i h 



90 
~dt 



h 2 A mD 2 
A6 + 



dfM 



2 m 



(q-q- 1 ) 



(42) 



where f q (p) = qp q 1 — q 1 p q 1 - 

(d) As a final example let us consider the K-entropy 



which can be obtained from Eq. (JJJ) after posing 

In K(p) = A ln {K} 



pl + K _ pi K 



2k 



dx 



(43) 



(44) 



The two quantities A = y/1 — k 2 and a = [(1 — k)/{1 + k)] 1 / 2k are constants and fulfil the 
relation (1 ± k) a ±K = A. 

The ^-logarithm is defined as 



2k 



(45) 



and reduces, in the k — ► limit, to the standard logarithm (ln {0} (x) = lnx). 
From Eq. we obtain the following continuity equation for the field p 



dp t-, 

9l + V ' 



D 



P 



— V (p 1+K + P 1 - K 



which coincides with that proposed in Ref. 0] 
The associated NSE is given by 

86 







(46) 



i h 



dt 



A6+ -mD 2 

2m Y 2 



+ A 



V In 



l 2 



1 + k) p K + (l-K)p- K 
+ V(aj)0, 



A In 



(47) 



and reduces to the Doebner-Golding Eq. (J32*j) in the k — »• limit as well as Eq. ({jH} reduces 
to the standard BG-entropy. 
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